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Abstract 


Using eight magnetic stations in low latitudes, hourly values 
of Dst are determined for the ICY. Three-hourly mean Dst values 
are compared with corresponding values based on Kertz's Dst indices. 
Three-hourly mean Dat is also compared with aj> indices. It is found 


that Dst correlates well with aj> in gross features and that during 
magnetic storms activity represented by a£ recovers more rapidly 
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Introduction 


Among transient variations of the earth's magnetic field the most 
outstanding in both magnitude and duration ,are those of magnetic storms. 
Though their form and severity vary considerably from one event to another, 
the majority of magnetic storms have certain common characteristics. 

One of the characteristics is the decrease of the horizontal force, 

H, all over the world in the equatorial and moderate latitudes. This 
decrease is very nearly axially symmetric and does not depend on longitude 
or local time. 

Superimposed on the general decrease in H there are variations in all 
three components that depend on local time; these variations are of large 
amplitude in the polar regions. 

It is thus convenient to analyze the storm variation D into two parts, 
Dst and DS. Given an instant of time Dst is the average of D over all 

longitudes. DS is defined by D - Dst. The two variations Dst and DS 

may be obtained for a specified latitude, or for a belt of finite width 
centered at some latitude. These definitions apply to all three magnetic 
components. 

For a well-defined storm Dst and DS can be determined as functions of 
storm time that is measured from the onset of the storm. 

Statistical analyses of magnetic storm variations have been made 
essentially in the manner described above by Moos (1910), Chapman (1918, 
1927, 1935, 1952), Vestine et al. (1947), Yokouchi (1957, 1958), Sugiura and 
Chapman (1960) and others. 
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However, the determination of Dst need not be limited to times of 
magnetic storms, but can be extended to periods of less magnetic activi'ty 
or even to magnetically quiet times. In fact, Dst may be determined 
continuously as a function of universal time regardless of occurrence of 
magnetic storms. Dst so determined should provide a measure of magnetic 
activity. 

Vestine et al. (1947) have published such continuous Dst; they gave 
hourly Dst curves for the Second Polar Year, September 1932 to August 1933 
(Carnegie Institution of Washington Publication 580 pp. 230-231). The 
hourly Dst values were obtained by averaging hourly departures from the 
annual mean for San Juan, Alibag and Honolulu, all in low geomagnetic 
latitude and spaced approximately evenly in longitude. 

More recently Dr. Akasofu and Professor Chapman determined Dst for 
several magnetic storms (c.f., Akasofu and Chapman, 1961), and their Dst 
curves have been extensively used by other workers, . Such Dst data are 
now in great demand by those who concern themselves with solar cosmic 
rays, solar plasma, earth's radiation belts and other related geophysical 
phenomena. 

Professor Chapman, Dr. Akasofu and the present author initiated 
determination of hourly Dst values for the IGY with a hope that such 
data will be published on a routine basis in the future. The work was 
carried out with assistance from Dr. J. C. Cain of NASA and Mr. E. J. 
Chernosky of Air Force Cambridge Research Laboratories. 

Dr. JCertz has derived three-hourly index for Dst for the IGY by a 
method different from the one adopted here. In Section 5, three-hourly 


Dst values obtained by our method are compared with the corresponding 
values based on Dr. Kertz's indices. 

2 • Magnetic Stations Used 

The selection of magnetic stations was made with the following 
considerations. (1) Stations in low latitudes should be used to 
eliminate as much as possible disturbances originating in the auroral 
zones. (2) Stations near the magnetic equator should be avoided, 
because the equatorial electrojet may introduce undesired variations; 
also magnetic storm variations are irregular and often augmented at the 
magnetic equator. 

TABLE 1 

Magnetic stations used - in the order of increasing 
geographic east longitude. 


Geomagnetic » Geographic 


latitude longitude ( E) latitude longitude ( E) 


Hermanus 

s 

33?3 

80?3 

S 

34° 

25' 

19° 

13 

Alibag 

N 

9?5 

143°7 

N 

18° 

38' 

72° 

53 

Kakioka 

N 

26?0 

206?0 

N 

36° 

14' 

140° 

11 

Apia 

S 

16?0 

26092 

S 

13° 

48' 

188° 

14 

Honolulu 

N 

21^0 

266°4 

N 

21° 

18' 

201° 

54 

San Juan 

N 

29°9 

3?2 

N 

18° 

23' 

293° 

53 

Pilar 

S 

20°3 

496 

S 

31° 

40' 

296° 

07 

M'Bour 

N 

2i 9 2 

55°1 

N 

14° 

24' 

343° 

03 
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(3) Stations should be distributed over longitude as uniformly as 
possible. (4) Reliable hourly H values are already available. 

The magnetic stations used are listed in Table 1 in the order of 
geographic longitude. Of the eight stations, five are in the northern 
and three in the southern hemisphere. The geomagnetic latitudes of the 
stations range from 9?5. to 33?3, 

3. Method of derivation of Dst 


For station i^, the observed magnetic field at universal time t 

may be considered as the sum of the permanent field H (including 
secular variation), solar daily variation , lunar daily variation 

and disturbance i.e. 


H (i) (t) = H 0 (i) (t) + Sq (1) (t) + L ( °(t) + D (l) (t) 


( L) 


Let the mean of H^^(t) over the 18 months from July 1957 to 

( 1 ) 

December 1958 be denoted by , and let this mean be a zero order 

approximation for H Let the deviation of HQ^*\t) from 

be A t) ; i.e. , 


(i) (i) d> 

H n (t) = H + &H (t) 

0 00 0 


( 2 ) 


We denote the deviation of observed H^*^(t) from H ^ ^ by H^^(t); 


i.e. , 


(i) (i) U) 

A H' (t) ■ H (t) - H qo 


(3) 


Substituting (1) and (2) in (3), we have 

AH (1) (t) = AH (1 \t) ♦ Sq (1) (t) ♦ L (1) (t) ♦ D (1) (t) (4) 
0 
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For hour t we take the mean of AH (i) (t) over the stations: 


AH(t) = £ H (t) + Sq( t) + L(t) + D(t) (5) 

where the upper bar in each term signifies the average over the stations. 
Since the lunar variation L is small,, we assume that the average 
of L is negligible, and ignore L altogether. Then we have 


A H(t) = a'h (t> + Sq( t) + D(t) (6) 

For station i_ for month M (where M = 1, 2, 18 for the 18 

months), we take the mean daily variation of the five international 


quiet days from Greenwich midnight to Greenwich midnight, and subtract 
from each of the 25 hourly mean values the linear change obtained by 
linearly connecting the two midnight values. This linear change is 
assumed to represent the non-cyclic change that should be subtracted 
from Sq. It should be noted that Sq is measured from midnight level 
instead of its daily mean level. 

The Sq variations so obtained were averaged over the stations. 

One would hope this mean Sq to be nearly zero, or at least very small. 

But our results showed that Sq was not comletely averaged out, and that 
the mean Sq had ranges of 8 to 19 gammas, a not quite negligible amount. 
The form of the residual Sq varied systematically with season. 

To represent Sq by a smoothly varying function throughout the whole 
period the mean Sq was expanded in a double Fourier series with month 
number M and Greenwich time T as two variables. 
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A™ cos (mT + cX ) cos ( nM+ 
n m 


f- 


) 


Coefficients up to n = m = 6 were determined. 

For each day mean Sq was synthesized with this series. In the 
computation a month number with one decimal was assigned to each day. 


The synthesized mean Sq was assumed to represent Sq ( t) in (6) , 


and was removed from AH(t), 


We are thus left with A IMt) + D ( t) . When this quantity 
4as plotted for the entire IGY, points remained near a fixed level 
during quiet periods and did not show any systematic drift all through 
the IGY. This was taken to mean that the secular variations were very 


nearly averaged out. Therefore, A H^(t) was assumed to be a constant. 
This constant determines the zero level for Dst* 

In determining the value of this constant it is important to 
exclude, as much as possible, days that are in the recovery phase of 
magnetic storms; inclusion of such days will lower the zero level of Dst. 
All the quiet periods in which two or more successive days had ap 


not exceeding 7 were selected, and the quantity A H (t) - Sq ( t) 


( = A Hq + D (t) was averaged over all these selected quiet days. 


This average, which was 31 gamma,, was taken to be the constant AH . 


Removing this constant from all the hourly values, we obtain D( t) . 
This average disturbance variation is the sum of Dst(t) and the average 


DS(t). Since the selected eight stations are all in low latitudes and 


are distributed approximately uniformly over longitudes, DS( t) is assumed 
to be negligible, and we take D( t) to represent Dst(t). 



7 


Summarizing, Dst is determined by 

Dst(t) = Hit) - H qo - Sq( t) - A H q (7) 

Here, the first term is the average over the stations of the hourly 

H value measured from tne mean H for the 1GY, the second term represents 

the mean Sq, and the third term is the mean Dst for the qufet periods. 

Since Dst varies with latitude, it is desirable to normalize 

to the variation appropriate to the equator. Assuming that the Dst 

field i s uniform and parallel to the geomagnetic axis, we multiply the 

hourly Dst values obtained in the manner described above by sec 8 , 

m 

where 0 is the mean geomagnetic latitude of the eight stations. In 
m 

our case 8 is 22° and the normalization factor is 1.08. 
m 

We will call the values normalized to the magnetic equator the 
equatorial Dst values. 

4, Tabulations and Graphs of Dst 

Hourly equatorial Dst values determined by the method described 
in the preceding section are tabulated in Table 2. In the tabulations, 
hourly values are given in units of gamma. 

The same equatorial hourly Dst values are graphically shown in 
Figure 1. Like Dr. Bartel's Kp diagram each line represents one solar 
rotation of 27 days. 

To compare our Dst with the three-hourly index ap, the hourly equatorial ^ 
Dst values were averaged over each of the same three-hourly intervals as 
for aj>, and the three-hourly Dst values so obtained are plotted together 
with aj> in Figure 2. Since Dst is negative during disturbed periods, the 
ap indices are plotted with negative sign. 
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Comparison of Dst with ag is made in Section 7. 

5. Comparison with Kertz's Dst Indices 

Dr. Kertz has determined Dst indices by a different method (see his 
article, pp. in this volume). 

To compare his results with ours, Dr. Kertz's indices were converted 
to three-hourly values in gammas. This conversion was made by multi- 
plying his indices by -3. 

In Figure 3, our three-hourly Dst (marked A ) and Dr. Kertz's Dst 
(marked B ) are plotted parallel to each other. 

Since Dr. Kertz took weighted overlapping means of three-hourly 
intervals, the curves representing his Dst values are on the whole smoother 
than ours. 

There appears to be a cftXirnal component not averaged out in Dr. 

Kertz's Dst. Examples are: September 10 to 20, 24 to 28, 1957; May 

22 to 24, 1958; October 14 to 17, 19 to 21, 1958; November 15 to 22 , 1958. 
These diurnal components are probably due to DS not being averaged out. 

The method we used does not involve any smoothing process. Hence 
our hourly Dst plots in Figure 1 show variations in greater detail than 
the plots in Figure 3 based on Dr. Kertz's results. For instance in 
Figure 1 an increase in H in the initial phase can be seen in many of 
the storms, 

6 . Pi scussions 

It is generally thought that the Dst field during a magnetic storm 
is due to a ring current. According to Chapman (1952) such a ring current 
was first suggested by Stormer in 1911, though Stbrmer's current ring was 
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not^/a complete ring, but a stream of electrons deflected round the earth 
on the afternoon side. Chapman further refers to A. Schmidt's study 
made in 1924 on an electrically neutral ring current. 

Chapman and Ferraro (1933, 1941) investigated a mechanism of form- 
ation of a ring current, its equilibrium, stability and decay. 

Singer (1957) proposed for the main phase decrease of the magnetic 

* 

field during a magnetic storm a mechanism, in which charged particles 
trapped in the earth's magnetic field drift longitudinally, thereby 
constituting net westward current. 

With stimulus from the discovery of the trapped radiation belts 
the problem drew renewed interest in the recent years. Thus Dessler 
and Parker (1959) gave a hydromagnetic interpretation of the effect of 
a ring current consisting of trapped particles. The magnetic field of 
a model ring current has been computed by Akasofu and Chapman (1961), 
and Akasofu, Cain and Chapman (1961). 

As to direct magnetic observation of a ring current by satellite 
or space probe there has been no definitive determination of the 
position of a ring current. However, Vanguard 3 magnetic measurements 
definitely indicate that the ring current must be above a few thousand 
kilometers from the earth’s surface. (Cain et al., 1962). 

We believe that the Dst field is essentially due to a ring current, 
and hence that our hourly Dst values represent the magnetic field 
variations of the ring current. The positive variations in Dst in the 
beginning of magnetic storms are due to the currents on the interface 
between the magnetosphere and solar plasma, but these are small compared 
with the ring current effect. 
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The hourly Dst plots as shown in Figure 1 demonstrate the 
development and decay of a magnetic storm ring current very clearly. 

Even when there is no well-defined magnetic storm there often appears 
to be weak ring current activity lasting continuously for a considerable 
length of time. Some of these variations may be due to polar disturb- 
ances not completely eliminated. 

The magnetic field observed on the earth's surface is modulated by 
the- strength of solar wind. Such an effect may be included in the Dst 
variations, 

7. Comparison of Dst with ap Indices 

Figure 2 gives a graphical comparison of three-hourly mean Dst 
with three-hourly magnetic indices ag. The latter indices are averages 
of three-hourly ranges for twelve stations lying between 48° and 63° 
geomagnetic latitude. The three-hourly Dst values presented here are 
based on hourly values of H for eight stations between 9°and 34° geo- 
magnetic latitude. 

The Dst field is due to the electric currents flowing at great 
distances from the earth, whereas the currents responsible for the magnetic 
disturbance represented by ag flow in the ionosphere. 

Though the two indices differ widely in their derivations and 
represent magnetic fields produced by currents flowing in so different 
altitudes, the curves in Figure 2 for these indices are remarkably parallel. 
Major Dst changes are accompanied by large ag variations of nearly equal 
magnitude. Correlation can often be seen between secondary peaks in Dst 
and those in ap . 



It is uncertain how far the details of our Dst are due to DS that 
is not completely eliminated. However, it is very unlikely that all the 
minor changes and irregularities in our Dst are due to DS not being averaged 
out. 

There is a tendency for ap to recover more rapidly than Dst after 
a magnetic storm. This means that in a magnetic storm, polar disturb- 
ance dies away more quickly than does the ring current. 

A Remark on the Induction in the Earth and in the Ionosphere 

If a magnetic field varying with time is applied to the earth, 
the magnetic induction within the conducting earth will modify the 
applied field. The effect is to increase the variation in H, and the 
magnitude of the effect depends on the time scale of the applied field. 

Since Dst involves variations with periods of wide range, the effect 
of the magnetic induction in the earth cannot be accurately estimated 
by a simple method. 

The earth, to some degree, is shielded by the ionosphere from varying 
magnetic fields produced by sources external to the ionosphere. This 
shielding effect depends on the rate of change of the applied magnetic 
field. 

To deduce the magnetic field of a ring current from our Dst values 
appropriate corrections must be made for the induced magnetic fields. 


Acknowledgements 


Professor S. Chapman initiated the publication of the hourly Dst 
data in the Annals of the IGY. I wish to express my thanks to him for 
his discussions and advice regarding the presentation ; to Dr. J. C. Cain 
for his advice on machine computation and for discussions; to Dr. S. -I. 
Akasofu for his participation in the original planning and discussions; 
and to Mr. E. J. Chernosky and Miss M. P. Hagan for their contribution 
in the laborious task of getting the observatory data On punch cards. 

I am greatly indebted to Mrs. Shirley Hendricks who performed the 
computation with an IBM 7090, and supervised plotting with an automatic 
plotter. 

Our sincere thanks are extended to the directors of the magnetic 
observatories who made the necessary data available to us. Assistance 
received from the World Data Center A for Geomagnetism is gratefully 
acknowledged* 



12 


References 

Akasofu, S. -I. , and S. Chapman 



1961 

J. Geophys. Research, 66, 1321 

Akasofu, S. 

•X . , J i C* 

Cain and S. Chapman 


1961 

J. Geophys. Research, 66, 4013 

Chapman, 8. 

19L8 

Proc. Roy. Soc . London, A, 95, 61 

«►- 


1927 

ibid. , 115, 242 


1935 

Terr. Mag. , 40, 349 


1952 

Annali di Geofisica, 5, 481 

Chapman, S. 

1952 

Ann. Geophys. , 8, 205 

Chapman S . , 

and V. C. 

A. Ferraro 


1933 

Terr. Mag. , 38, 79 


1941 

ibid. , 46 , l 

Dessler, A. 

J . , and E. 

N. Parker 


1959 

J. Geophys. Research , 64, 2239 

Cain, J. C. 

, 1. R. Shapiro, J. D. Stolarik and J. P. Heppner 


1962 

J. Geophys. Research, 67, 5055 

Moos, N. A. 

F. 1910 

Magnetic Observations Made at the Government 
Observatory, Bombay, for the Period 1846 to 
1905, and their Discussion, Part II. The 

Phenomenon and its Discussion, Bombay. 

Singer , S. 

F. 1957 

Trans. Am. Geophys. Union, 38, 175 


Sugiura, M. , and S. Chapman 

1960 Abh. Akad. Wiss. Gottingen, Mathr Phys. 

Klasse, Sonderheft Nr, 4, Gttttingen 

Vestine, E. H. , L. Laporte , 1, Lange and W, E. Scott 

1959 The Geomagnetic Field, Its Description and Analysis, 
Carnegie Institution of Washington Publication 580 , 
Washington, D.C. 

Yokouchi, Y. 1957 Mem. Kakioka Observatory, 8, (1), 27 
1958 ibid . ,8, (2) , 13 








DATS IN SOLAR ROTATION INTERVAL 



Figure I coni. : HOURLY Dst PLOTS 
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Figure I cont. : HOURLY Dst PLOTS 
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